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I.  INTRODUCTION 


Passage  of  a  shock  moving  with  uniform  velocity  over  a  thin  layer  of 
heated  fluid  results  in  the  formation  of  a  thermal  precursor  (Fig.  1).  This 
interaction  was  studied  analytically  in  Ref.  1  and  experimentally,  for  weak 
shocks,  in  Refs.  2  and  3.  More  recently,  an  inviscid  numerical  code  was  used 
to  evaluate  thermal  precursor  flow  fields.  An  analytical  model  developed  in 
Ref.  5  yields  results  that  are  in  approximate  agreement  with  numerical 
examples  reported  in  Ref.  4. 

The  present  study  provides  an  estimate  for  the  turbulent  boundary  layer 
development  within  the  precursor  region.  Results  from  Refs.  4  and  5  are  used 
as  a  guide  to  define  the  inviscid  flow  external  to  the  boundary  layer.  The 
turbulent  boundary  layer  within  the  precursor  is  then  evaluated  using  an 
unsteady  momentum  integral  formulation.  Boundary  layer  thickness  parameters 
are  defined  in  Appendix  A.  Related  problems,  namely,  the  turbulent  boundary 
layer  induced  by  passage  of  a  shock  wave  over  a  semi-infinite  plate,  the 
turbulent  wall  boundary  layer  induced  In  a  shock  tube,  and  the  turbulent 
boundary  layer  induced  by  the  Impulsive  motion  of  a  semi-infinite  plate,  are 
discussed  in  Appendices  B  and  C. 
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II.  THEORY 


The  thermal  precursor  flow  field  is  described.  A  solution  is  then 
presented  for  the  turbulent  boundary  layer  induced  within  the  precursor. 

A.  FLOW  DESCRIPTION 

The  inviscid  thermal  precursor  flow  field,  illustrated  in  Fig.  1,  occurs 
when  the  stagnation  pressure  downstream  of  a  shock  of  strength  =  us/a4 
region  4  is  less  than  the  static  pressure  downstream  of  the  incident  shock 
Mg  =  Ug /a j  In  region  1.  In  the  latter  case,  a  portion  of  the  fluid  in  region 
2  moves  forward  and  drives  a  shock  in  the  thermal  layer  with  velocity 
u  >  u_.  The  interface  between  the  fluids  from  regions  1  and  4  is  indicated 

S  L  S 

by  a  dashed  line  in  Fig.  1.  The  quantities  *Sp,  x^,  and  xgt  in  Fig.  1  denote, 
respectively,  the  location  of  a  local  wall  stagnation  point,  the  interface 
location  at  the  wall,  and  the  location  of  the  shock  wave  within  the  thermal 
layer.  The  incident  shock  location  is  denoted  xg.  The  velocities  corre¬ 
sponding  to  xspJ  Xi,  and  xgt  are  denoted  ugp,  u^,  and  ugt,  respectively.  The 
results  of  Ref.  4  indicate  that,  after  an  initial  transient,  these  velocities 
tend  to  constant  values. 

The  boundary  layer  within  the  precursor  region  is  illustrated  in  Fig.  2. 
The  dashed  line  again  denotes  the  interface  between  fluids  from  regions  1  and 
4.  The  boundary  layer  in  the  intervals  x^  <  x  <  xgt  and  x  <  0  can  be  reduced 
to  a  steady  state  by  considering  coordinate  systems  wherein  the  thermal  layer 
shock  and  the  free  stream  shock  are  stationary,  respectively.  The  latter  are 
termed  “shock  induced  boundary  layers"  in  Fig.  2  and  can  be  accurately 
evaluated  (e.g..  Ref.  6).  The  boundary  layer  in  the  region  0  <  x  <  Xj 
cannot  be  reduced  to  a  steady  state.  The  solution  is  complicated  by  the 
presence  of  the  interface  within  the  boundary  layer.  A  solution  for  these 
regions  is  given  in  the  next  section.  A  similar  boundary  layer  is  associated 
with  the  impulsive  motion  of  a  semi-infinite  plate,  with  a  shock  tube  wall, 
and  with  shock  passage  over  a  semi-infinite  plate. ^  References  7  to  9 
consider  laminar  boundary  layers.  The  corresponding  unsteady  turbulent 
boundary  layer  problems  are  discussed  in  Appendix  B. 
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SOLUTION  OF  BOUNDARY  LAYER  EQUATIONS 


The  present  flow  model  Is  Indicated  In  Fig.  3.  It  is  assumed  that  each 
of  the  steamwise  stations  xg,  xgp,  x^,  and  xgt  moves  with  constant  velocity. 
The  separation  between  xg  and  xgp  is  ignored  (Fig.  3a).  The  latter  assumption 
is  used  to  simplify  the  solution  for  the  initial  portion  of  region  III  and 
does  not  affect  regions  I  and  II.  The  flow  field,  at  time  t  *  t',  is 
illustrated  in  Fig.  3b.  The  free  stream  velocity  in  regions  I  and  II  is 
assumed  to  equal  the  interface  velocity  u^ ,  in  approximate  accord  with  the 
results  of  Refs.  4  and  5.  The  streamwise  variation  of  fluid  properties  of 
state,  external  to  the  boundary  layer  in  regions  I  and  II,  is  ignored.  The 
effect  of  the  latter  assumption  is  discussed  in  "Concluding  Remarks." 

The  momentum  integral  equation  for  each  of  the  regions  I,  II,  and  III,  in 
Fig.  3,  can  be  expressed  in  the  form*® 
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where  6**  and  8  are  as  defined  in  Appendix  A.  Wall  shear  is  related  to  the 
local  boundary  layer  thickness  6  by  the  Blasius  relation^ 
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Here,  state  properties  are  based  on  a  mean  temperature 


0.5  (£u  l)  +  0.22  (3l-  1)  (3) 


where  Tw  and  Tf  denote  wall  and  recovery  temperatures,  respectively.  (See 
Appendix  A.) 
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We  further  assume  that  6  /6  and  8,.  do  not  vary  with  x  or  t  in  each  of 

the  regions  cf  interest. 

Equations  (1)  and  (2)  can  then  be  combined  to  yield 


The  boundary  layer  thickness  6  has  been  assumed  to  be  zero  at 
(Fig.  3b).  The  boundary  conditions  are  then 
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f(£gt)  =  0  (7b 

Since  streamwise  variations  of  A  and  are  being  neglected,  Eq.  (6a)  can  be 
integrated  to  yield 


f  -  A  +  C  a  -  5*) 


(8 


where  C  is  a  constant. 

The  boundary  conditions  [Eqs.  (7a)  and  (7b)]  provide  two  separate  values 

for  C.  The  resulting  two  solutions  for  f  can  be  joined  at  station  £  such 
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that  f,  and  therefore  6  are  continuous  at  £  .  The  result  is 
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where 
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Incompressible  flow 


(9c 


The  ability  to  satisfy  the  boundary  conditions  at  £gp  and  £gt  and  match 

ft  ft 

at  5  -  £  is  a  consequence  of  the  fact  that  £  is  a  singular  point  of  Eq. 

m 

(6a).  The  boundary  layer  thickness  is  a  maximum  at  (Fig.  3).  The  surface 
shear  is  a  minimum  at  that  station. 

The  boundary  layer  thickness  and  local  surface  shear  coefficient  can  be 
expressed  in  physical  variables.  The  result,  for  £gp  <  €  <  Cm»  is 
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and  for  5  <  €  <  € 
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where 


♦  -  (pm/pJA/5  <w1/5 


The  ratio  9/6  in  Eqs.  (10)  and  (11)  has  been  normalized  by  the  value  for  incom¬ 
pressible  flow.  The  solution,  however,  is  not  limited  to  incompressible  flow. 


In  the  limit  £  “0,  Eqs.  (10a)  and  (10b)  correspond  to  the  steady-state 

sp 

solution  for  a  turbulent  boundary  layer  on  a  semi-infinite  flat  plate. ^  When 

£  *  0,  Eqs.  (10a)  and  (10b)  account  for  streamwise  motion  of  the  boundary 

sp 

condition  6  ■  0  at  xgp.  Equations  (11a)  and  (lib)  correspond  to  the  solution 
for  a  turbulent  boundary  layer  behind  a  shock  moving  with  uniform  velocity. ^ 
The  solution  is  accurate,  in  region  I,  if  the  appropriate  free  stream 
conditions  are  used. 


Region  III 


The  boundai /  layer  in  region  III  corresponds  to  the  boundary  layer  behind 

the  incident  shock  which  moves  with  velocity  u„  ■  x_/t.  Let  subscript  2 

s  s 

denote  free  stream  conditions  behind  the  incident  shock  and  define  £  *  x/u2t. 

A  derivation,  similar  to  that  for  Eqs.  (11a)  and  (lib)  yields  for  £  <  5 
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where  4  -  (p  /p,)4^5  (v  /v.)175,  £*  »  (0/6**)  and 
m  4  m  z  z  L 
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Here,  ^27pi  *s  t*ie  density  ratio  across  the  incident  shock.  Equations 
(12a)-(c)  are  applicable  in  the  region  x  <  0  and,  for  the  case  ug  -  ugp  in  the 
region  0  <  x  <  x 

sp. 


Time  Dependence  at  Fixed  Station 


Equations  (10)  through  (12)  can  be  put  into  a  form  which  is  convenient 
for  evaluating  boundary  layer  properties  at  a  fixed  station  (e.g.,  x'  in  Fig. 
3a).  Assume  that  the  interaction  starts  at  time  t  -  0  (Fig.  3a)  and  let  tg, 


.  <  «r.  r.  * 


tSp,  tif  and  tgt  denote  interface  arrival  times  at  x'.  Boundary  layer 
properties  at  station  x'  are  found  from  Eqs.  (10)  through  (12)  by  use  of  the 
following  substitutions: 


_1_ 

u 

e 


x  - 


uspc 

177 


(13a) 


“at*  -  x 


(t  -  St^i 
Cm  -  Cst 


(13b) 


1  u  t  -  x  t  -  t 
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where  tm  -  tj/£  .  It  is  clear  from  Eqs.  (13a)  and  (13b)  that  the  boundary 
layer  solution  is  continuous  at  time  t  *  tffl. 

C.  LAMINAR  BOUNDARY  LAYER  CASE 

The  present  procedure  can  be  used  to  treat  the  case  of  a  laminar  boundary 
layer  induced  by  a  thermal  precursor  and  related  laminar  problems In 
particular,  the  quantity  6^*  in  Eqs.  (4a)  and  (5b)  is  replaced  by  62  and  the 
constant  A  is  expressed  as 


A  “  2  x  S/(p  u  2) 
w  e  e 


(14) 
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2 

The  solution  for  f  *  6  /t  is  then  given  by  Eq.  (9).  The  maximum  boundary 
layer  thickness  occurs  at 
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IS 


** 


11, 


(15) 


0.384  (incompressible  ) 


which  may  be  compared  with  the  value  £  ■  7/9  for  an  incompressible  turbulent 

m 

boundary  layer. 


III.  CONCLUDING  REMARKS 


A  solution  has  been  presented  for  the  turbulent  boundary  layer  within  a 
thermal  precursor.  The  solution  is  limited  because  differences  in  fluid  state 
properties  between  regions  I  and  II  have  not  been  taken  into  account.  An 
improved  theory  can  be  developed  by  the  piecewise  application  of  Eq.  (8)  in 
regions  xgp  <  x  <  xm,  xm  <  x  <  x^,  and  x^  <  x  <  xgt  wherein  the  boundary 

condition  5  -  0  is  satisfied  at  xgp  and  xgt  and  6  is  continuous  at  xffl  and 

x^.  This  approach  does  not  appear  warranted  without  an  improved  model  for  the 
boundary  layer  in  region  xm  <  x  <  x^. 

The  numerical  thermal  precursor  solutions  of  Ref.  4  indicate  that,  after 
an  initial  transient,  u.  =  ugt«  Thus  the  separation  distance  xgt  -  x^  remains 
constant.  The  fixed  separation  xgt  -  Xj  is  a  consequence  of  the  fact  that  the 
fluid  which  crosses  the  shock  in  the  thermal  layer  is  not  constrained  to 
remain  between  xg£  and  x^  (as  in  a  shock  tube),  but  is  ejected  from  the  wall 
region  near  x^.  The  flow  is  then  similar  to  supersonic  flow  (of  the  heated 

layer  gas)  over  a  blunt  body  (cold  precursor  gas).5  Thus,  after  long  times 

xg(;  -  x^  is  small  compared  with  xgt  -  xgp.  Hence  region  I  is  small  compared 
with  region  II.  In  this  case  (i.e.,  £  ♦  £  ),  the  turbulent  boundary  layer 
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in  the  precursor  can  be  estimated  from  Eqs.  (10)  and  (11)  by  using  free-stream 
properties  from  region  II  In  each  of  these  equations. 

The  boundary  layer  behind  a  shock  moving  over  a  serai-infinite  plate,  the 
shock  tube  wall  boundary  layer,  and  the  boundary  layer  associated  with 
impulsive  fluid  motion  over  a  serai-infinite  plate  are  illustrated  in  Figs.  4 
through  6,  respectively.  These  flows  represent  special  cases  of  Eqs.  (10)  and 
(11)  and  are  discussed  in  Appendix  B.  Numerical  results  for  the  shock-induced 
turbulent  boundary  layer  on  a  semi-infinite  plate  are  given  in  Ref.  12  and 
compared  with  the  present  model  in  Appendix  C.  Reference  12  presents  surface 
shear  and  boundary  layer  thickness  for  the  case  of  standard  air  with  Mg  “  1.58 
and  xg  ■  419.3  cm.  These  results,  based  on  a  numerical  mixing  length  model, 
are  given  in  Fig.  7.  For  this  case,  the  present  model  predicts  that  the 
maximum  boundary  layer  thickness  and  the  minimum  surface  shear  occur 


Flow  field  associated  with  passage  of  shock  over 
semi-infinite  flat  plate.  Dashed  line  denotes  particles 
that  were  at  x  -  0  at  time  t  ■  0 


Fig. 


.  Wall  boundary  layer  in  shock  tube  for  case  of  thii 
expansion  wave  and  ideal  (unperturbed)  centerline 
Dashed  line  denotes  particles  at  diaphragm  locatli 
(  x  »  0)  at  time  t  *  0 


Flow  field  associated  with  impulsive  start  of  fluid 
motion  over  a  semi-infinite  flat  plate.  Dashed  line 
denotes  particles  that  were  at  x  »  0  at  time  t  «*  0 


u  — 

CO  II 


O  t;  w* 

c 

O)  CO  * 
GO  CM 

«c  • 

<j  e  « 

4J 

^  to  GO 
O  iji 

U-(  ^ 


0.150 

0.149 

0.148 

0.147 

0.146 

0.145 

0.144 

0.143 

CD  CD  CD  CD  O 

CO  Cvi  n — ■  CD) 

JUJDI'IQ  g/L(0^/0n|  Uq'L 

UI3 

UV3HS  BOVdanS  1V301 

SS3N)l3IHi  H3AV1  AyV0N(108 

APPENDIX  A 


CHARACTERISTIC  BOUNDARY  LAYER  THICKNESSES 


Let  subscript  e  denote  local  free  stream  properties  in  a  wall  fixed 
coordinate  system.  Boundary  layer  defect  thicknesses  are  then 


<5  /<5  =  /  [l  -  (pu/ogue)]dy/6 
o 

1 

9/<$  “  /  (Pu/p  u  )[  1  -  (u/u  )]dy/6 
e  e  e 

o 

**  1 

<5  /6  *  /  (p/pg)[l  -  (u/ug)]dy/<5 

o 


where  6  and  6  are  displacement  and  momentum  thicknesses,  respectively.  The 
quantity  6  appears  in  Eq.  (1)  and  reduces  to  6*  for  the  case  of  an 
incompressible  fluid. 

Equations  (A-la)-(c)  can  be  evaluated  for  a  turbulent  boundary  layer  if 
velocity  and  density  profiles  are  known.  A  suitable  choice  is^ 


u/ug  «  (y/6) 


P„  T. 


+  -c(^] 


u  '  '  u  ' 

e  e 


where 


b  =  (T  /T  )  -  1 
r  w 


c  =  [<T  /T  )  -  l]  (T  /T  ) 
w  r  e  e  w 


T  /T  -  1  +  [(y  -  l)/2]r(o)M2 
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APPENDIX  B 


RELATED  BOUNDARY  LAYERS 


The  turbulent  boundary  layers  Induced  by  passage  of  a  shock  over  a  semi- 
infinite  flat  plate,  by  shock  tube  flow  and  by  the  impulsive  start  of  fluid 
motion  over  a  semi-infinite  flat  plate,  have  features  that  are  similar  to  the 
boundary  layer  in  a  thermal  precursor.  These  boundary  layers  are  now 
discussed . 

The  flow  associated  with  shock  passage  over  a  semi-infinite  plate  is 

illustrated  in  Fig.  4.  The  dashed  line  indicates  particles  that  were  located 

at  x  *  0  at  time  t  *  0.  A  solution  for  the  laminar  boundary  layer  case,  in 

the  limit  Mg  -*■  1,  is  presented  in  Ref.  9,  in  which  unsteady  two-dimensional 

laminar  boundary  layer  equations  were  considered.  The  flow  in  Fig.  4  is 

similar  to  the  flow  in  Fig.  3b  with  u„„  “  0.  Hence,  a  momentum  integral 

sp 

solution  for  the  turbulent  boundary  layer  case  leads  to  Eqs.  (10)  and  (11) 

with  £  =  u  *  0  in  Eqs.  (10a)  and  (10b).  The  flow  in  region  0  <  x  <  x_ 

sp  sp  111 

corresponds  to  a  steady  flow  over  a  semi-infinite  flat  plate.  The  flow  in 
region  xm  <  x  <  xg  corresponds  to  the  boundary  layer  induced  by  a  moving  shock 
wave.  Station  x  (i.e.,  C  *  6/5  )  represents  the  station  at  which  6  is 
matched  for  these  two  flows  and  is  a  maximum.  In  this  case,  the  assumption  of 
uniform  free-stream  properties,  external  to  the  boundary  layer,  is  exactly 
satisfied. 

The  wall  boundary  layer  induced  in  a  shock  tube  is  illustrated  in  Fig. 

5.  A  finite  difference  solution  for  the  boundary  layer  in  region  0  <  x  <  x^ 
was  presented  in  Ref.  8.  for  the  case  of  laminar  flow.  Equations  (10)  and 
(11)  can  be  used  to  estimate  the  turbulent  boundary  layer  everywhere  in  the 
shock  tube  if  the  expansion  wave  thickness  is  neglected  (as  illustrated  in 
Fig.  5)  and  if  fluid  property  differences  between  the  driver  gas  and  the 
driven  gas  are  neglected.  The  boundary  layer  in  region  x  <  x  <  x  is  found 
from  Eq.  (10)  with  ugp  and  ?gp  replaced  by  ugx  and  5gx  =  x^/x^,  respec- 


tively.  Note  that  u_„  and  £  are  both  negative.  The  boundary  layer  in 

ex  ex 

region  x^  <  x  <  x  is  obtained  from  Eq.  (11).  It  is  seen  that  the  maximum 
®  *  * 
value  of  6  again  occurs  at  £  ■  8/6 

m 

The  flow  associated  with  impulsive  motion  of  flow,  over  a  semi-infinite 

flat  plate,  is  illustrated  in  Fig.  6.  The  flow  is  assumed  to  be  impulsively 

started  at  time  t  =•  0.  The  dashed  line  represents  particles  that  were  at 

x  =  0  at  time  t  *  0.  The  laminar  case,  discussed  in  Ref.  7,  contains  a 

momentum  integral  solution  and  a  discussion  of  the  singular  nature  of  the  full 

laminar  boundary  layer  equations.  A  momentum  integral  solution  for  the 

turbulent  case  yields  Eqs.  (10a)  and  (10b)  with  £  *  u  “0,  for  0  <  x  <  x_. 

sp  sp  U1 

Hence,  here  too,  the  flow  is  equivalent  to  steady  flow  over  a  semi-infinite 

plate.  For  x  >  xffl,  the  boundary  layer  thickness  and  shear  are  time  dependent 

but  do  not  vary  with  x.  Boundary  layer  properties  are  found  by  evaluating 

Eqs.  (10a)  and  (10b)  with  £  -  u  «  0,  at  x  *  x_  *  (6/5  )(u  t).  The 

sp  sp  m  e 

boundary  layer  for  x  >  xm  is  the  turbulent  analog  for  the  laminar  Rayleigh 
problem^  wherein  an  infinite  plate  is  set  into  motion  at  time  t  ■  0.  The 
present  procedure,  to  match  at  x^  is  similar  to  the  approach  used  in  Ref.  7 
for  the  laminar  case. 

In  all  of  these  cases,  as  in  the  case  of  the  thermal  precursor,  the 

**. 

maximum  value  of  6  and  the  minimum  value  of  wall  shear  occur  at  £  0/6 

m 

As  previously  noted,  the  latter  represents  a  singular  point  which  joins  the 

similarity  solution  for  £  <  £  to  the  similarity  solution  for  £  >  £  and.  is  in 

m  m 

accord  with  the  theory  of  Stewartson.^ 


APPENDIX  C 


COMPARISON  WITH  REFERENCE  12 


Numerical  turbulent  boundary  layer  results  are  reported  in  Ref.  12  which 
provide  the  variation  of  surface  shear  and  boundary  layer  thickness  behind  a 
shock  moving  over  a  serai-infinite  plate  (Fig.  4).  A  mixing  length  model  was 
used  and  finite  difference  calculations  were  performed.  The  results  of 
Ref.  12  are  compared  herein  with  the  present  model. 

The  local  shear  and  boundary  layer  thickness  variation,  from  Ref.  12,  is 
given  in  Fig.  7.  These  data  are  for  the  case  of  air  and  Ms  -  1.58,  Tj  »  290 
K,  and  P^  ■*  1  atmosphere.  The  data  represent  values  at  the  instant  of  time 
when  the  shock  is  at  xg  «  419.3  cm.  For  this  case  6*/6  -  0.1330,  0/6  «  0.09916, 
6**/6  *  0.12845,  £  -  2.0,  £,  -  1.0,  and  £  i  0/6**  -  0.772.  In  the  vicinity 

of  £  *  £  ,  the  numerical  results  for  6  and  Cf  in  Ref.  12  are  approximately  equal 
to  1.8  and  0.8  times  the  values  deduced  herein,  respectively. 

The  station  at  which  surface  shear  is  a  minimum  and  boundary  layer 

thickness  is  a  maximum  is  of  primary  interest.  The  present  model  predicts  this 

location  to  be  at  £  *  0.772,  which  is  included  In  Fig.  7.  It  is  seen  that  this 

in 

prediction  is  in  close  agreement  with  the  numerical  results  of  Ref.  12.  This 
agreement  tends  to  confirm  the  present  model. 
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LABORATORY  OPERATIONS 
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The  Aerospace  Corporation  functions  as  an  "architect -engineer"  for 
national  security  projects,  specializing  In  advanced  military  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  invest Igat ions  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  is  the  technical  staff's  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 


J.N 

£ 


y. 

?• 


Aerophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamlnat ion,  thermal  and  structural 
control;  high  temperature  thermoraechanlcs ,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  excimer  laser  development  including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out -of -f ield-of -view  rejection, 
applied  laser  spectroscopy ,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  faul t -tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications; 
microwave  sem t conductor  devices,  microwave/millimeter  wave  measurements , 
diagnostics  and  radiometry,  microwave/millimeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysts  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy-induced 
env  i  moments . 

Space  Sciences  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomy, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
i nst  r umentat ion. 
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